Model of a VT symmetric Bose-Einstein condensate in a delta-functions double well 
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The observation of VT symmetry in a coupled optical wave guide system that involves a complex 
refractive index has been demonstrated impressively in the experiment by Riiter el al. (Nat. Phys. 
6, 192, 2010). This is, however, only an optical analogue of a quantum system, and it would be 
highly desirable to observe the manifestation of VT symmetry and the resulting properties also 
in a real, experimentally accessible, quantum system. Following a suggestion by Klaiman et al. 
(Phys. Rev. Lett. 101, 080402, 2008), we investigate a VT symmetric arrangement of a Bose- 
Einstein condensate in a double well potential, where in one well cold atoms are injected while in 
the other particles are extracted from the condensate. We investigate, in particular, the effects of 
the nonlinearity in the Gross-Pitaevskii equation on the VT properties of the condensate. To study 
these effects we analyze a simple one-dimensional model system in which the condensate is placed 
into two VT symmetric (^-function traps. The analysis will serve as a useful guide for studies of the 
behaviour of Bose-Einstein condensates in realistic VT symmetric double wells. 

PACS numbers: 03.65.Ge, 03.75.Hh, 11.30.Er 



Beginning with the seminal paper by Bender and 
Boettcher in 1999 [1], parity-time (VT) symmetric quan- 
tum mechanics has attracted ever increasing attention 
over the past decade because it offers a class of complex 
Hamiltonians which, in spite of their non-Hermiticity, 
possess discrete real energy eigenvalue spectra. More- 
over, these Hamiltonians feature the property of branch 
points, i.e., the coalescence of both energy values and 
eigenfunctions when some parameter in the Hamiltonian 
is varied, a phenomenon impossible in Hermitian quan- 
tum mechanics (but known to appear for resonances in 
the continuous spectrum, see, e.g., [2J]). 
Recently VT symmetry has been realized experimentally 
in structured optical waveguides [3, |4j , where the com- 
plex index of refraction is manipulated by introducing 
loss and gain terms. These experiments make use of the 
quantum-optical analogy that the wave equation for the 
transverse electric field mode is formally equivalent to the 
one-dimensional Schrodinger equation. It would, how- 
ever, be desirable to observe VT symmetry also in a real 
quantum system. 

Klaiman et al. [5| have suggested a quantum scenario 
analogous to the waveguide experiments in which a Bose- 
Einstein condensate is placed in a double well potential, 
and loss and gain is realized by removing atoms in one 
well and coherently adding particles in the other. These 
authors pointed out that to have a close analogy with 
the optics experiments the nonlinearity in the Gross- 
Pitaevskii equation governing these condensates should 
be kept small. Here we want to ask the opposite question: 
What are the effects of the nonlinearity on the VT sym- 
metry on such an arrangement of a Bose-Einstein conden- 
sate? It is namely exactly the nonlinearity, proportional 
to ^(x)! 2 , in the Gross-Pitaevskii equation that compli- 
cates matters. A necessary condition for the Hamilto- 
nian to be VT symmetric is that the imaginary part of 
the potential is an odd function, and the real part an 



even function of x. The latter cannot be assumed from 
the outset for \yb(x)\ 2 when solving the Gross-Piatevskii 
equation. 

In this paper we will investigate the effects of the non- 
linearity on the VT symmetry in the spirit of a model 
calculation by considering the situation where the dou- 
ble well is idealized by two delta-function traps, with 
loss added in one trap and gain in the other. We will 
demonstrate that the stationary solutions of the Gross- 
Pitaevskii equation indeed preserve the VT symmetry of 
the nonlinear Hamiltonian, and merge in a branch point 
at some critical value of the loss and gain, beyond which 
the symmetry is broken. Our results show that it will 
be a worthwhile enterprise to investigate VT symmetric 
Bose-Einstein condensates in realistic double well poten- 
tials, and possibly pin down physical parameters where 
VT breaking could be observed in a real experiment. 
A model which mimics the physical situation of a BEC in 
a symmetric double well with loss and gain has already 
been investigated by Graefe et al. 6-8] in the frame- 
work of a two-mode Bose-Hubbard-type VT symmetric 
Hamiltonian. As an optical analogue, in the two-mode 
approximation Ramezani et al. [9[ have recently looked 
at a mathematical model of a VT symmetric coupled 
dual waveguide arrangement with Kerr nonlinearity. It is 
one objective of this paper to see which features of these 
models are recovered when actually solving the nonlinear 
VT symmetric Gross-Pitaevskii equation. 
For a system where a real delta-function potential is aug- 
mented by a VT symmetric pair of delta-functions with 
imaginary coefficients, bound states and scattering wave 
functions have been calculated by Jones [10( . His interest 
was devoted to the quasi- Hermitian analysis of the prob- 
lem, and no nonlinearity was present. Jakubsky and Zno- 
jil [ll| have considered the explicitly solvable model of a 
particle exposed to two imaginary VT delta-function po- 
tentials in an infinitely high square well, and determined 



the energy spectrum. The nonlinear Schrodinger equa- 
tion for a delta-functions comb was studied by Witthaut 
et al. [12J with the aim of gaining insight into the proper- 
ties of nonlinear stationary states of periodic potentials. 
Also, there exists a vast amount of literature on soli- 
tons and Bose-Einstein condensates in periodic optical 
and nonlinear lattices with VT symmetry and their non- 
linear optical analogues (see, e.g., 13l-l2l|). But to the 
best of our knowledge the basic problem of two VT sym- 
metric delta-function double wells with Gross-Pitaevskii 
nonlinearity has not been considered so far. 
The Gross-Pitaevskii equation we analyze in this paper 
has the form 
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with n G C, Re(«;) > 0, and 7 real. 

It consists of two delta-function traps with distance a, 
located at b = ±a/2, with a real attractive part of the 
potential and imaginary gain/loss terms whose strengths 
are determined by the parameter 7, and a nonlinear term 
with amplitude g, which arises from the contact interac- 
tion of the condensate atoms. Units have been chosen 
in such a way that the strength of the real part of the 
delta-function potential is normalized to unity. While 
the delta-function potentials are VT symmetric, it is not 
clear a priori that the equation itself is VT symmetric 
since this requires the nonlinear term to be a symmetric 
function. 

For vanishing nonlinearity, we find that the simple quan- 
tum mechanics model captures, for both eigenvalues and 
wave functions, all the effects of a VT symmetric wave 
guide configuration in optics. This is essentially due to 
the fact that two attractive delta-function potentials have 
exactly two bound states which correspond to the two su- 
permodes in the wave guide arrangement. 
For nonvanishing nonlinearity we have solved the Gross- 
Pitaevskii equation ([I]) numerically using a procedure 
in which the energy eigenvalues are found by a five- 
dimensional numerical root search. The free parameters 
which have to be adjusted in such a way that a physically 
meaningful wave function is obtained are the eigenvalue 
k as well as initial conditions for the wave function and 
its derivative. Since the overall phase is arbitrary we can 
choose it such that \l/(0) is a real number. Therefore five 
real parameters remain, viz. the real part of ^(0), and 
the real and imaginary parts of both ^'(0) and k. Phys- 
ically relevant wave functions must be square integrable 
and normalized. The normalization is important since 
the Gross-Pitaevskii equation is nonlinear and the norm 
influences the Hamiltonian. This gives in total five condi- 
tions which have to be fulfilled: The real and imaginary 
parts of \1/ must vanish for x — > ±00, and the norm of 
the wave function must fulfill ||^|| — 1 = 0. 
Outside the delta-function traps the Gross-Pitaevskii 
equation ([!]) coincides with the free nonlinear 
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FIG. 1. (Color online) Eigenvalues k, of the nonlinear equa- 
tion {TJ as functions of the size of the loss/gain parameter 7 
for a — 2.2. The value of the nonlinearity parameter chosen is 
g — 1. The case g — is drawn (dashed lines) for comparison. 
A branch of complex conjugate eigenvalues appears which bi- 
furcates from the ground state branch before the critical value 
of 7 where the branches of the real eigenvalues coaelesce in an 
exceptional point. There a pair of complex eigenvalues only 
emerges after an analytical continuation of the nonlinearity 
in the Gross-Pitaevskii equation. 



Schrodinger equation, which has well known real 
solutions in terms of Jacobi elliptic functions (cf., e.g., 
121 |22L |23[). The function which solves the equation 
in the ranges \x\ > b for the attractive nonlinearity 
considered here and decays to zero for \x\ — > 00 is 
ch(kx, 1) = l/cosh(Kx). We find that once the cor- 
rect eigenvalues and eigenfunctions are obtained our 
numerical wave functions exactly show this behaviour. 
Fig. Q] shows the results for the eigenvalues k calculated 
for a value of the nonlinearity parameter g — 1.0 and a 
trap distance of a = 2.2 as functions of 7. The results 
for the case g = are also shown for comparison. It can 
be seen that even with nonlinearity there still exist two 
branches of real eigenvalues up to a critical value j CI ss 
0.4, at which the two eigenvalues coincide. There also 
appears a branch of two complex conjugate eigenvalues, 
but surprisingly these are born, not at 7 cr , but at the 
smaller value of 7 « 0.31 where they bifurcate from the 
real eigenvalue branch of the ground state. This implies 
that there is a range of 7 values where two real and two 
complex eigenvalues coexist. 

At this point it is useful to establish a link with the model 
of a VT symmetric Bose-Hubbard dimer with loss and 
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FIG. 2. (Color online) Real and imaginary parts and moduli 
of the eigenstates of the nonlinear Hamiltonian in equation 
(1) for g — 1, o = 2.2, and 7 = 0.35, (a) ground state and (b) 
excited state. The wave functions are VT symmetric, and the 
moduli are symmetric functions, producing the VT symmetry 
of the total nonlinear Hamiltonian. 



FIG. 3. (Color online) Real and imaginary parts and moduli 
of the eigenstates of the nonlinear Hamiltonian in equation (1) 
for g = 1, a = 2.2, and 7 = 0.5, (a) solution with imaginary 
part of k > and (b) imaginary part of k < . The VT 
symmetry is broken, the moduli are not symmetric functions, 
and the VT symmetry also of the nonlinear Hamiltonian is 
broken. 



gain investigated by Graefe et al. [8|. An eigenenergy 
spectrum with a structure similar to the one in Fig. [T] 
also appeared in their calculations (see Fig. 13 in [8|). In 
the model, stationary states correspond to fixed points 
of the motion of a vector on the surface of the Bloch 
sphere, whose types can be classified according to the 
eigenvalues of the Jacobian matrix. In the region where 
only two real eigenvalues exist the solutions correspond 
to centers, while in the region with four eigenvalues the 
solutions correspond to a center and a saddle point, and 
a sink and a source. The center and saddle point collide 
at the branch point and vanish. This behaviour is in 
complete agreement with the results shown in Fig. [TJ It 
may be concluded that the familiar branching scheme 
known for VT symmetric Hamiltonians quite generally 
will be changed into a scheme of the type as shown in 
Fig. [T] if a nonlinearity is added to the Hamiltonian. 

Fig. [21 shows the real and imaginary parts of the ground 
state and the excited state determined numerically for 



g = 1, a = 2.2, and 7 = 0.35, below the the critical value 
7 cr ~ 0.4. The VT symmetry of each wave function 
is evident since their real parts are even functions and 
their imaginary parts odd functions of x. From the VT 
symmetry of the wave function follows that the modulus, 
also shown in Fig. [5] is an even function, and with it the 
nonlinear term in equation (fT]). We therefore have the 
important result that the nonlinear Hamiltonian picks 
as eigenfunctions exactly those states in Hilbert space 
which render the nonlinear Hamiltonian VT symmetric! 
In the ground state, which emerges from the symmetric 
real wave function for 7 = g = 0, the symmetric contri- 
bution from the real part still dominates, while for the 
excited state, which originates from the antisymmetric 
solution for 7 — g = 0, the antisymmetric contribution 
from the imaginary part prevails. 

The VT symmetry of the wave functions is broken for 
the eigenstates with complex eigenvalues. Fig. [3] shows 



as an example the wave functions obtained for g = 1, a = 
2.2, and 7 = 0.5 for the corresponding pair of complex 
conjugate eigenvalues re. It can be seen that the real 
and imaginary parts are no longer even or odd functions, 
and therefore VT symmetry is lost. As a consequence, 
the moduli of the wave functions also are no longer even 
functions of x. Thus we find that beyond the branch 
point not only the VT symmetry of the wave functions 
is broken but also that of the nonlinear Hamiltonian! 
For the states with complex eigenvalues there is, how- 
ever, an important difference between the case with and 
without nonlinearity: for complex eigenvalues the modu- 
lus squared of the wave functions grows or decays propor- 
tional to exp(— 2Im(re 2 )£), and so does the nonlinear term 
in equation (1). Therefore strictly speaking the solutions 
presented here can only describe the onset of the tempo- 
ral evolution of the two modes, for times Im(re 2 )i <C 1. 
To obtain the full time evolution in the region beyond the 
branch point one would have to solve the time-dependent 
Gross-Pitaevskii equation. This clearly goes beyond the 
scope of our simple model calculation. In Fig.[3]the mode 
with positive imaginary part of re is the one which begins 
to decay, as expected it is more strongly localized in the 
trap with loss, while the mode with negative imaginary 
part is the one which starts to grow and is more strongly 
localized in the trap with gain. 

The fact that at the branch point two real solutions co- 
alesce without giving rise to two solutions with complex 
eigenvalues contradicts the usual behaviour seen at ex- 
ceptional points. Obviously these solutions cannot be 
found by solving the nonlinear Gross-Pitaevskii equation 
in its form ([T]), but require an analytical continuation of 
the nonlinear Hamiltonian beyond the critical point j cr . 
The reason is that the nonlinear term g\^>\ 2 is a non- 
analytic function, and some care has to be taken when 
analytically continuing the Hamilitonian beyond the ex- 
ceptional point. 

In the VT symmetric regime up to 7 cr we have ^*(x) — 
\f(— x). Therefore on the way to the bifurcation point we 
can write the nonlinearity for the VT symmetric states 
in the form ^^(x)! 2 = g' , S(x)'i>(~ x). This function can 
be continued analytically. In the numerical calculation 
the additional condition j ^>(x)^>(— x)dx = 1 must be 
introduced to fix the phase of the nonlinearity in the 
VT broken regime. In the (then) six-dimensional root 
search also Im('I'(0)) must varied. As a result we find two 
more complex conjugate solutions that emerge from the 
coaelescing states, see Fig. [TJ These states are not VT 
symmetric, and no longer possess vanishing imaginary 
parts at the origin. 

In this paper we have analyzed the simple quantum me- 
chanical model of a Bosc-Einstein condensate in VT sym- 
metric delta-function double traps by directly solving the 
nonlinear Gross-Pitaevskii equation. We find two sta- 
tionary eigenstates with real eigenvalues which at a crit- 
ical value of the loss/gain parameter merge in a branch 



point. We have the important result that the wave func- 
tions are VT symmetric. As a consequence their moduli 
are even functions, and therefore the nonlinear Hamilto- 
nian selects as solutions exactly such states which make 
itself VT symmetric. We also find a branch of two com- 
plex conjugate eigenvalues for which the VT symmetry 
of the wave functions is broken, and with it that of the 
nonlinear Hamiltonian. 

An unexpected result is that, with the nonlinearity 
present, the branches of complex conjugate eigenvalues 
do not bifurcate from the point where the real eigen- 
values conincide, but emerge at a smaller value of the 
gain/loss term from the eigenvalue branch of the ground 
state. On the other hand, at the critical value of the 
gain/loss parameter we find the behaviour characteristic 
of a branch point, i.e. the coalescence of both eigenval- 
ues and eigenfunctions, but no pair of complex conjugate 
eigenvalues seems to emerge. These are found only after 
continuing analytically the nonlinear term in the Gross- 
Pitaevskii equation. Note, however, that as stated before, 
for complex eigenvalues the squared modulus of the wave 
function becomes time-dependent, and a description us- 
ing the stationary Gross-Pitaevskii equation breaks down 
anyway. This does, however, not affect the main result of 
our paper, namely the existence of VT symmetric eigen- 
functions and the VT symmetry of the Hamiltonian also 
when the nonlinearity is present. 

We have considered the case of an attractive nonlinearity 
but found that the same behaviour occurs for repulsive 
nonlinearity. 

The results of our model calculation make one to expect 
that similar VT behaviour should also prevail in Bosc- 
Einstein condensates in more realistic double wells [24( 
with VT symmetry, in one or more dimensions. Also 
in addition to the nonlinearity resulting from the short- 
range contact interaction condensates with a long-range 
dipole-dipole interaction [25| could be considered. In- 
vestigations of the Gross-Pitaevskii equation in these di- 
rections are under way. It would also be interestin g to 
extend the quasi-Hermitian analysis given by Jones [10( 
and to investigate whether for the nonlinear VT sym- 
metric Hamiltonian considered in the present paper the 
construction of a metric operator is possible with respect 
to which the nonlinear Hamiltonian is quasi-Hermitian. 
Furthermore it would be worthwhile looking for simple 
matrix models which show the behaviour of the eigenval- 
ues found for finite nonlinearity. 

We thank Eva-Maria Graefe and Miloslav Znojil for help- 
ful comments. 
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